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We study a model where two scalar fields, that are subdominant during inflation, decay into 
radiation some time after inflation has ended but before primordial nucleosynthesis. Perturbations 
of these two curvaton fields can be responsible for the primordial curvature perturbation. We 
write down the full non-linear equations that relate the primordial perturbation to the curvaton 
perturbations on large scales, calculate the power spectrum of the primordial perturbation, and 
finally go to second order to find the non-linearity parameter, /nl- We find large positive values 
of /nl if the energy densities of the curvatons are sub-dominant when they decay, as in the single 
curvaton case. But we also find a large /nl even if the curvatons dominate the total energy density 
in the case when the inhomogeneous radiation produced by the first curvaton decay is diluted by 
^ . the decay of a second nearly homogeneous curvaton. The minimum value min(/NL) = —5/4 which 

^ ' we find is the same as in the single-curvaton case. 



3 ■ 



> 
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I. INTRODUCTION 



Theories beyond the standard model often contain a large number of scalar fields in addition to the standard-model 
fields. In the very early universe it is natural to expect the initial values of these fields to be displaced from the 
minimum of their potential. If they are displaced by more than the Planck scale then they can drive a period of 
inflation. But if they are displaced from their minimum by less than the Planck scale they will oscillate about the 
minimum of their potential once the Hubble rate, H, drops below their effective mass. An oscillating massive field has 
the equation of state (averaged over several oscillations) of a pressureless fluid. Thus the energy density of a weakly 
i— — i, interacting massive field tends to grow relative to radiation in the early universe. Such fields must therefore decay 
before the primordial nucleosynthesis era to avoid spoiling the standard, successful hot big bang model. And if the 
energy density of a late-decaying scalar is non-negligible when it decays then any inhomogeneity in its energy density 
will be transfered to the primordial radiation p| [2] • This is the curvaton scenario for the origin of structure [!, 0, Hf • 
A curvaton field, \, is supposed to have a negligible energy density during inflation but once the Hubble rate drops 
£SJ . below the curvaton mass after inflation, the curvaton energy density grows relative to radiation, reaching its maximum 
value, f2 Xi decay, just before the curvaton decays. If all the species are in thermal equilibrium and the baryon asymmetry 
is generated after the curvaton decays, then curvaton mechanism generates adiabatic density perturbations. If not 
then the curvaton can leave a residual isocurvature perturbation @, 0| correlated with the curvature perturbation [f| . 
\ The amplitude of isocurvature modes are severely constrained by current data [H, [13, [0, G3 • Thus if, for example, 
■ the baryon asymmetry is produced by the out-of-equilibrium curvaton decay, then we require fi Xl decay ~ 1- 

If we take seriously the multiplicity of scalar fields in the early universe then we should consider models where 
more than one field can contribute to the primordial density perturbation on large scales (l3| . Several authors have 
considered the combination of perturbations from a curvaton field and the inflaton field driving inflation [lil [Tol [i"6| . 
More recently Choi and Gong (l7j considered the primordial perturbations that may result from multiple curvaton 
fields, showing that the presence of more than one curvaton field affects the amplitude of residual isocurvature 
perturbations and their correlation with the curvature perturbation. In this paper we extend their analysis to study 
the non-linear curvature perturbation and how the multiple late-decaying scalar fields may affect the non-Gaussianity 
of the primordial curvature perturbation. For simplicity we restrict our analysis to the case of two curvatons, but it 
should be straightforward to extend our analysis to three or more curvatons. 

Deviations from an exactly Gaussian distribution of the primordial density perturbation is conventionally given in 
terms of a non-linearity parameter, /nl [HI- The current upper bound from the WMAP three- year data is |/nl| < 
114 while Planck is expected to bring this down to |/nl| < 5 [HI]. Galaxy cluster surveys can offer complementary 
constraints [20| . Measurement of /nl would give a valuable test of inflation. If primordial perturbations originate 
from fluctuations in a canonical inflaton field, driving slow-roll inflation, then /nl is less than unity [HI, [23. However, 
if primordial perturbations originate from fluctuations in a single curvaton field, then /nl ~ decay [1, US HH HE 
Hy, [2?| • Therefore /nl could thus be large if the curvaton does not dominate the energy density of the universe when 
it decays. In this paper we shall show that when we consider the decay of two curvaton fields it is also possible that 
/nl is large even when the energy densities of both curvaton fields are dominant when they decay. 

This paper is organized as follows. In Sec. [IT] we introduce a non-linear definition of curvature perturbations, define 
the perturbation power spectrum and bispectrum, and the non-linearity parameter /nl which describes the non- 
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Gaussianity of perturbations at leading order. In Sec. IIIII we write down the full non-linear equations that relate the 
curvaton perturbations to the total and radiation perturbations at the first-curvaton decay and at the second-curvaton 
decay, and finally to the primordial curvature perturbation. Then we solve these equations up to the second order 
and find the resulting /nl of the primordial perturbation before nucleosynthesis. The general expression for /nl is a 
non-trivial function of four parameters. To get an insight into its behavior we analyze some special cases in Sec. IIVI 
Finally, in Sec. [V] we present concluding remarks. As the notation becomes rather heavy in the two-curvaton scenario, 
we list most of the symbols used in this paper in Appendix [AJ 



II. PRELIMINARIES 



A. Primordial curvature perturbation 



The primordial density perturbation can be described in terms of the non-linearly perturbed expansion on uniform- 
density hypersurfaces [H[ (see also [29l l3(il] ) 

1 M*> x ) do 

C(t,x) = ^,x) + - / (1) 
3Jp(t) P+P 

where N = J Hdt is the integrated local expansion, p the local density and p the local pressure, and p is the 
homogeneous density in the background model. 

We will expand the curvature perturbation at each order (n) as 

C(*> x ) = Z!^C(n)(*,x), (2) 

n=l 

where we assume that the first-order perturbation, is Gaussian as it is proportional to the initial Gaussian field 
perturbations. Higher-order terms describe the non-Gaussianity of the full non-linear £. 

Working in terms of the Fourier transform of £, we define the primordial power spectrum as 

(C(k 1 )C(k 2 )) = (2vr) 3 P c (fc 1 ) ( 5 3 (k 1 + k 2 ) . (3) 

The average power per logarithmic interval in Fourier space is given by 

47rfc 3 
(2tt) 3 

and is roughly independent of wavenumber k. 
The primordial bispectrum is given by 



W = ^(W- (4) 



(C(k 1 )C(k 2 )C(k 3 )} = (27r) 3 B(k 1 , k 2 , fc 3 )5 3 (k 1 + k 2 + k 3 ) . (5) 

The bispectrum vanishes for a purely Gaussian distribution, and hence is non-zero only at fourth and higher-order. The 
amplitude of the bispectrum relative to the power spectrum is commonly parameterized in terms of the non-linearity 
parameter, /nl, defined such that [l8[ 

B(k 1 ,k 2 ,k 3 ) = (6/5)/NL[P(fci)P(fc 2 ) + 2perms] . (6) 

Higher order statistics, like trispectrum (see e.g. [13, Hi], H2j) or mn non-linear probability density function of the 
primordial £ (27j . can give also valuable information on non-Gaussianity in the curvaton model, but in this paper we 
consider the bispectrum only. 



B. Curvaton perturbations 

We will consider the primordial curvature perturbation produced by the decay of two scalar fields a and b. Without 
loss of generality, we assume that the curvaton a decays first when H = T a followed by the decay of the curvaton b 
when H = Th, where Tj, <T a . 

Small-scale (sub-Hubble) vacuum fluctuations of any light scalar field are stretched by the expansion to large (super- 
Hubble) scales during a period of inflation in the early universe. If the curvaton fields are weakly-coupled, massive 
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scalar fields whose masses are less than the Hubble scale, 3> m, during inflation, then they acquire an almost 
scale-invariant spectrum of perturbations on super-Hubble scales, 



H * x 



V Sa , = V 6b . = [-£) , (7) 



where is the Hubble rate at Hubble exit. In Rcfs. 33, 34, 35j, |36j, |37| it has been shown that in general, if slow-roll 
conditions are satisfied, the non-Gaussianity of field perturbations at Hubble exit is small (at least for the three-point 
and four-point correlators). Consequently, in what follows, we assume that the field perturbations at Hubble exit 
have Gaussian and independent distributions, consistent with weakly-coupled isocurvature fluctuations of light scalar 
fields [Mill. 

Once the Hubble rate drops below the mass of each curvaton field the field begins to oscillate. The local value 
of each curvaton can evolve between Hubblc-cxit during inflation and the beginning of the field oscillations. We 
parameterize this evolution by functions g a and gb, but we assume the two fields remain decoupled so that their 
perturbations remain uncorrelated. Thus at the beginning of curvaton oscillations, the curvaton fields have values 

a sc = 9a (a*) (8) 
bosc = 9b(K) ■ (9) 

We can define non-linear perturbation for each curvaton analogous to the total perturbation (|TJ) 

( a = SN+hn(^)=SN+hn(^) . (10) 
3 \Pa J 3 \a osc J 

where we take the energy density to be proportional to the square of the field value, p a oc a^ sc , when the field begins to 
oscillate, and similarly for (b- On uniform-curvaton-density hypersurfaces (40] | £ a and are curvature perturbations. 
The curvature perturbation £ Q (and £b) becomes constant on scales larger than the Hubble scale once each curvaton 
starts oscillating, and while we can neglect energy transfer to the radiation, T a < H. Thus we require 

T a < m a and r b < m b , (11) 

consistent with the assumption that we are dealing with weakly coupled fields. For simplicity we will further assume 
that r a <C mb and that the energy density of both curvatons is negligible when they begin oscillating, which ensures 
that we can neglect any gravitational coupling between fields in the overdamped regime before they begin oscillating. 

At first order the curvature perturbation (| 10[) for each curvaton can thus be given in terms of the field perturbations 
on spatially flat hypersurfaces (6N — 0) 

, 2 5a 2g> 

Ca(1) = H^^"" (12) 
= =if^ (13) 

3 b osc 3 9b 

where we have Taylor expanded the functions g a and gb, defining g' a = dg a /da*, and g' b = dgb/db*. 

The power spectra P^ a and P^ a(1) are the same at leading order. The higher order corrections are generally so 
small that we use P^ a and f£ o(1) interchangeably in the following. Therefore the power spectra of £ a and £fc when the 
curvatons start to oscillate are related by 

Pa - P 2 p<« . (14) 

where 

0=4^- as) 

9 a/ 9a 

If we assume linear evolution between Hubble exit and the beginning of curvaton oscillations, this factor reduces to 
the ratio of the background curvaton field values at Hubble exit, 

0* 
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At second order Eq. (fit))) gives on spatially flat hypersurfaces (SN = 0) 

3 L 9a9 

2 ^ g'i 



Co(2) = ~\ ( 1 " ) Ca(l) . ( 17 ) 



Cb(2) - -| f 1 - ^pr) Cfc(i) , (18) 

where £ a (i) and Cbfi) are given in ([12]) and tfl3]). Analogous to ([6]) we can define the non-linearity parameters of each 
curvaton and /^ L . Then comparing to (fl"7]) and (fT8)) we find 

5 SaSa 

4 I ff« 



/nl = -tU-^M - (19) 



/nl = -\[^- 9J f)- ( 2 °) 

If the evolution of field values is linear between Hubble exit and the beginning of curvaton oscillations, we have 
/nl = /nl = ~5/4- Although it is possible to construct potentials which lead to a non-linear evolution [4l[, we 
assume henceforth in this paper that = = 0, consistent with weakly interacting fields. 



C. Calculating /nl 

Before calculating how the two curvaton perturbations contribute to the primordial curvature perturbation in detail, 
we can first consider the general form of the primordial power spectrum and bispectrum that will result from a generic 
model where the local values of two Gaussian fields determine the primordial perturbation. 

The primordial curvature perturbation after both curvatons decay can be written up to second order in terms of 
the first order curvaton perturbations as 

C = C2 = C2(l) + ^2(2) 

= A(a{l) + BCb(l) + |CC a (l) + |-DCb(l) + 2-ECa(l)C&(l), (21) 

where £2 is the total curvature perturbation after the second decay. The first order part £ 2 (i) is Gaussian, since it 
is a linear function of the field perturbations at Hubble exit (and for weakly coupled light scalar fields these field 
perturbations are Gaussian). The second order part C2(2) is non-Gaussian, since it contains squares of Gaussian 
variables. A, B, C, D and E are coefficients that will depend on background quantities at the time of the first and 
second decay. In Sect. HTT1 we derive the expressions for them in the two-curvaton decay scenario. 
The primordial power spectrum at leading order is 



P ( = A 2 P Ca + B 2 P Cb . (22) 



Employing the relation (|14[) we can write 



= ^^t- (23) 

= WTpm p <' (24) 

Since the two curvatons are uncorrelated with one another, i.e., 

<Ca(i)(ki)C K i)(k2))=0 (25) 

for any wave vectors ki and k 2 , we find that at leading order the three-point correlator of the primordial perturbation 
is 



C(k 1 )C(k 2 )C(k 3 )) = 

_±A 2 c(Ca(l)(kl)Ca(l)(k 2 )^ 

+|ASs/c a( i)(k 1 )C 6 (i)(k 2 )(C a( i)*C6(i))(k 3 ))l + 2 permutations of {k l7 k 2 ,k 3 }, (26) 
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where * denotes a convolution. For example, 



(Ca(i)*Ca(i))(k 3 ) = / ^qCa(i)(q)Ca(i)(k3 -q). (27) 



Note that there are two ways to form two pairs of correlators in the first two terms on the right-hand-side of (|26p . 
Therefore, the three-point correlator simplifies to 

(C(k 1 )C(k 2 )C(k 3 )) = {[A 2 CP Cam (h)P Ca(1) (k 2 ) + B 2 DP Qb{1) {k 1 )P Qb{1) {k 2 ) 

+\ABEP Ul) {k 1 )P Cb(1) {k 2 )} + 2 perms} x (2n) 3 5^ (kx + k 2 + k 3 ) . (28) 

Substituting and JM]) into (J2HD leads to 

C (ki C k 2 C k 3 - ^ 

{P c (fc 1 )P c (fc 2 )+ 2perms} x (27r) 3 ^ 3 )(k 1 +k 2 + k 3 ). (29) 
Hence the bispectrum, defined in Eq. ([5]) , is in the form given by Eq. ([6]) , where the non-linearity parameter is 



5 CA 2 + \(5 2 EAB + [3 A DB 2 

6 (A 2 + P 2 B 2 ) 2 



/nl = ~ a /' . -r-rd • (30) 



III. FULL NON-LINEAR EQUATIONS 

We will estimate the primordial density perturbation produced by the decay of two curvaton fields some time after 
inflation has ended using the sudden-decay approximation Q , generalizing the non- linear analysis of Ref. [13] to the 
case of two curvatons. In this approximation the curvaton fields and the radiation are treated as non-interacting fluids 
except at the instant of decay when all the energy density of the curvaton is transfered to radiation. This should be 
a good approximation for scales much larger than the decay time, r , and has been shown to give a good estimate 
of the primordial non-Gaussianity when compared against numerical simulations [HI, [H, [HJ • 

Before either curvaton has decayed, the radiation (from inflaton decay products) has a definite equation of state, 
p-y — (l/3)p 7 , and so do the oscillating curvaton fields, p a = Pb = 0. Thus we have three non-interacting fluids with 
barotropic equations of state and hence three curvature perturbations |l} which are constant on large scales [io| 

C 7 - 5N + \ In , (31) 



4 Vp 7 



3 \Pa 

Cb = SN + \ In I . (33) 



3 \Pb 

On the spatial hypersurface where H — T ai there is an abrupt jump in the overall equation of state due to the 
sudden decay of the first curvaton into radiation, but the total energy density is continuous 

P 7l l + Pbl = P70I + Pal + Pbl ■ (34) 

Here p 7l i is the radiation energy density immediately after the first curvaton decay, pbi is density of the second 
curvaton at the time of the first curvaton decay, p l0 \ is the density of pre-existing radiation just before the first 
curvaton decay, and p a \ is the density of the first curvaton just before the first decay when it is converted to radiation. 
Eq. ([34]) simplifies to 

P71I = P7ol + Pal ■ (35) 



Note that the first-decay hypersurface is a uniform-density hypersurface and thus, from Eq. ([T]) the perturbed 
expansion on this hypersurface is SN — £1, where Ci denotes the total curvature perturbation at the first-decay 
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hypersurface. Hence on the first-decay hypersurface we have 



p-rol 




(36) 


Pill 


= p 7ll e 4 ^), 


(37) 


Pal 


= Paie 3 ^-^, 


(38) 


Pbl 


= p bl e 3 «»-^ , 


(39) 



where we employed Eqs. (f3TH3"3"|) . In the curvaton scenario we usually assume the initial radiation is homogeneous 
before the curvaton decays, so that £ 7o = 0, but in this section we keep £ 7o m our calculations for generality. On 
a uniform-density hypersurface the total energy density is homogeneous. Therefore an infinitesimal time before the 
first decay we have 

P~f„l + Pal + Pbi = Pi- (40) 
Substituting here Eqs. (|36[) . (f38|) and ([39]) . and dividing by the total energy density p\ we end up with 

tt 7o ie 4(c -o- Cl) + n al e 3 ^-^ + n 61 e 3 K»-W = 1 , (41) 

where f2 7o i = p 7o i/pi, = Pal/ Pit and O&i = Pbl/ Pi are the energy density parameters of the pre-existing 
radiation, the first curvaton, and the second curvaton at the first curvaton decay, respectively. On the other hand, an 
infinitesimal time after the first curvaton has decayed into radiation the total energy density can be written as 

Aril +Pbi = Pi , (42) 

which gives an equation 

7ll e 4 ^i- Cl ) + fW ( ^- Cl) = 1 . (43) 

The total curvature perturbation £ is the perturbation on the decay hypersurface and thus is continuous between the 
two phases. However the energy density of radiation changes abruptly at the decay time and the radiation curvature 
perturbation is discontinuous, (~ (0 7^ £ 7l . 

Similarly, an infinitesimal time before the second decay we have 

fi 7l2 e 4(C ^ ~ C2) + ft fc2 e 3(C ^ C2) = 1 • (44) 

while an infinitesimal time after the second decay we have 

fi 722 e 4(c ^- C2) - 1. (45) 

Above f2 7l 2 and £l 2 are the radiation and the second curvaton energy density parameters just before the second 
curvaton decay, f2 72 2 is the radiation density parameter just after the second curvaton decay, and £2 is the total 
curvature perturbation at the second-curvaton decay hypersurface. As the radiation is the only constituent left after 
the second decay, the radiation curvature perturbation £ 72 equals the the total curvature perturbation £2, and Eq. 1(15]) 
reduces to ri 72 2 = 1 as we will see in the next subsection. After both curvatons have decayed we have a single radiation 
fluid with equation of state p — p/3 and hence the curvature perturbation remains constant. Thus we identify the 
primordial curvature perturbation, 

C = C 72 - (46) 

We can Taylor expand the exponential functions, e x = 1 + x + x 2 /2 + . . ., in Eqs. (|1T1 - |4"5")) to express the final 
curvature perturbation, £, in terms of the initial curvaton perturbations, Ca and at any given order. 

A. Zeroth order 

Taylor expanding Eq. (|4ip at zeroth order gives 

^ 70 i + n»i + flu = 1 . (47) 
This is identically true, since we are studying spatially flat models with fi tot = 1. Similarly Eq. (|43|) gives 

n 7ll + n bl = i. (48) 
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Therefore fi 71 i is a redundant variable, and in what follows we can replace it by 1 — Qbi- To zeroth order, equation 
© is 

tt 7l2 + fi b2 = 1 • (49) 

Therefore f2 7l 2 is a redundant variable, and in what follows we can replace it by 1 — fl 2- Finally Eq. (f4"5"l) is identically 
true as it states f2 72 2 = 1. At all higher orders Eq. (|45|) states C 72 — C2- After the second decay only radiation is left. 



B. First order 



1. First decay 



At first order, Eq. (|4Tj) gives 

4^7oi [C 70 (i) -Ci(i)]+3fi i[C„ ( i) -Ci(x)]+3Jl6i[Ci(i) -Ci(i)] =0- (50) 
From this we solve for the total curvature perturbation, Ci(i)i a ^ the first decay surface. It is 

Cl(l) = / 7 olC 7o (l) + /olCo(l) + fblCb(l) ) (51) 

where 

/701 = 4O 70l + 3fili + 3fi 6 i ' (52) 

/. 3r2 n i . . 

Ial ~ 4fi 70l + 3ft Ql + 3fi M ' 1 } 

p _ 3f^i . . 

Jbl 4fi 7ol + 3n al + 3fl bl ' 1 ' 

These first order curvature transfer efficiency parameters at the first decay obey the relation / 7o i + f a \ + f^i = 1. In 
what follows we use this to eliminate / 7o i- 

After the first curvaton decays, but before the second curvaton decays, the curvature perturbation in the radiation, 
C 7l , will remain constant on large scales. Eq. (|43|) relates the radiation curvature perturbation and the total curvature 
perturbation immediately after the first decay. At first order it reads 

4fi 7ll [c 7l( i) - + 3n M [& ( i) - Cm)] - 0, (55) 

which gives 

C 7 i(i) = -RiCi(i) - (Ri - 1) Ch(i) j (56) 

where 

4-Obi 

Rl = 4-4^ (5?) 

3 + fal (58) 



3(1 - Ai) + fal 



Note that if the density of the second curvaton is negligible when the first curvaton decays then we have R\ = 1 and 

C 7 i(i) = Ci(i)- 

Finally, we substitute ([51]) into (JSJ 



C 7l (l) = Rl (1 - fal - fbl) C 70 (l) + RlfalCad) + [1 - J2l (1 - /bl)] Cfc(l) ■ (59) 

This remains constant between the first and second decay and hence, at the second decay we can use this as the 
incoming radiation perturbation. 
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2. Second decay 



Following the same procedure as above, we find from Eq. (|44[) at first order (or equivalently using (|5ip but dropping 
curvaton 6, then relabeling a — > 6, 1 — > 2 and — > 1) 

C2(l) = /7i2C 7 l(l) + fb2(b(l) , (60) 

where 

^ = 40 7 ^+ 3fi 62 ' (61) 

/b2 = 3 ^ b2 j-g 2 ) 
4J7 7l 2 + 3f2b2 

These first order curvature transfer efficiency parameters at second decay obey the relation / 7l 2 + fa = 1, and hence 
Eq. (HOD reads 

C2(l) = (1 - hi) C T i(l) + /fc2Cfc(l) ■ (63) 

Finally using Eq. for C7i(i) we obtain 

C2(l) = Rl{l~ fal - /bl) (1 - fa) Cro(l) + r aCa(l) + ?"bCb(l) 1 ( 64 ) 



where 



Rlfal (1 - Aa) (65) 
(l-/«)(3 + /«l)/„l 



m-M+u ' (66) 

r h = 1 - R t (1 - fa) (1 - fa) (67) 

_ (1 - fbl)fa(3 + fal) + fblfal /r> q \ 

- 3(f-/ M ) + / al • (68) 



From Eqs. (|43|) and (|4^|) we have £2(1) = Ct 2 (i) = C(i)- Therefore, if the pre-existing radiation perturbation vanishes, 
£70(1) = ^> ^ e P r i morc iial first order curvature perturbation after the second decay is 

C(i) = r Co(i) + nCb(i) ■ (69) 

Then in Eq. (|2f p at first order we can identify 

A = r Q , and B = r b , (70) 
and the power spectrum for the primordial curvature perturbation is given by Eq. 



^=rl + frl. (71) 

Ca(l) 

where /3 2 defined in Eq. (|f 4p gives the ratio between the initial power in the curvaton a and that in the curvaton b. 

Note that our result differs from that presented recently by Choi and Gong [I?) due to the presence of the extra 
factor Ri which arises due to the difference between the uniform total density hypersurface and the uniform radiation 
density hypersurface when curvaton a decays, if the density of the curvaton b is not negligible. When = we find 
from Eq. (|58j) that R\ = 1 and we recover the simpler result [TtJ 

C72U) = (1 _ /h2)/olCo(l) + fb2(,b(l) ■ (72) 

The transfer coefficients r a and r& in this case are shown by the thick solid lines in Figs. [Hand El 

However in general when 7^ and J?i ^ 1 we find a novel effect where the inhomogeneous density of the 
curvaton b may lead to a perturbation in the radiation density after the curvaton a decays, even if the curvaton a 
is homogeneous, see Eq. (15^1) . A local overdensity of the curvaton b delays the decay of curvaton a due to a local 
gravitational time dilation, resulting in a local overdensity in the radiation after decay. In practice fa > (3/4) fa so 
this is usually a small correction. The dotted lines in Figs. [T] and [2] show this effect for the case fa = fa/2- 
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FIG. 2: Same as Fig.[T]but now for log 10 rt- The contours of equal r;, are shown, from left to right, for rt = 10 3 , 10 2 ' 5 , 10 2 , 
lGT 1 - 6 , 10" 1 , and 10" ' 5 . 
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When considering r a and r& it should be noted that the range of fa, fa, and fa is [0, 1], but the allowed range 
of fa is constrained for given / al and fa. Since f a i + fa = 1 — / 70 i, and / 7ol 6 [0, 1], we always have a trivial 
constraint 

fa < 1 - /al ■ (73) 

Furthermore, we can rewrite /f,i in a form 



4 - fi a i — Ofti 4 - f2 



61 

Since 4 3i Q b is an increasing function of time, and fa — 4 3 _^ , we find 

The multiplier of / 2 simplifies to 1 + / a i/3. So we have a constraint 

/m<(1 + /«i/3)/ w . (76) 

For example, when fa — fa/2 the constraint (|73")l forbids a small region in the upper right corner of Figs. [1] and [5] 
where fa ~ 1 and fa - 1. 

Equality in Eq. (|T6[) corresponds to the case where f2f,i = f2 2, i.e., the second curvaton decays at the same moment 
as the first curvaton. Then the transfer efficiencies simplify to r a = fa, and r& = (1 + /ai/3)/ 2 = fa, indicated by 
the thin solid lines in Figs. [1] and [2 (Note therefore that the case where fa = fa describes a situation where the 
curvaton b decays some time after the curvaton a.) 

C. Quasi-second order 

In the next section we present the full second order calculation of £, but here we consider a simplified treatment, 
which will give the correct form of /nl generated by multiple curvaton decays in the limit where the non-Gaussianity 
is large. 

In this quasi-non-linear approximation we will use the linearized form for the full curvature perturbation (|32J) 
and (|33[) in terms of the density perturbations on spatially flat hypersurfaces: 

t = (77) 

3 Pa 

C, = ^, (78) 
3 pb 

but we will include non-linear contributions to the curvaton densities. The energy density of the curvaton a is 
p a oc a 2 = (a + da) 2 and hence we have 

Pa a \a J 

For simplicity we neglect any non-linear evolution of the curvaton field [i.e., we take = g' b ' = in Eq. ([5]) and ©J 
so that the curvaton field perturbations can be taken to be Gaussian and we have 

3 



4'«(i) 

where 



Ca = Cad) + jQl) . (80) 



r 2Sa * (sn\ 

Ca « " 3^7' (81) 



and similarly for the curvaton b. 
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Now we use the first order result for the primordial curvature perturbation produced by the two curvaton decays 
but instead of Ca(i) an( i Cfc(i) we use the quasi-nonlinear £ a and Cb of Eqs. (|77[) and Q78[) . Then we have from Eq 



C = ^aCa + r b Cb 



(1) 



Ca(l) 



,6(1) 



(82) 



Hence, comparing this result with Eq. (|2ip. we have in our quasi-second-order approximation 



A = 


r a , 


(83) 


B = 


n, 


(84) 


C = 


3 

2 ra ' 


(85) 


D = 


3 

7 b ' 


(86) 


E = 


0. 


(87) 



Substituting these into Eq. (|30[) the nonlinearity parameter reads 



/•quasi 

'NL — 



4^3 
b 



(3 4 r 



We derived this result using only first-order formulas for £, but including the full second order expression for the 
energy density. The true second order Q includes corrections of order £ 2 , but we expect the above result to give a 
good approximation when £( 2 ) ^ Cfi) t i- e -j when the non-Gaussianity is large. In the next section we present the full 
second-order calculation. 



D. Second order 

1. First decay 

At second order, Eq. ([4T|) reads 

tt 70 i[4(c 7o( i) - Cm))] 2 + n 0l [3(c tt(1) - Cm))] 2 + n w [3(C6 ( i) - Cm))] 2 

+4fi 7ol [c 7o(2) - Ci( 2) ] + 3r> al [c a (2) - Ci( 2) ] + 3n bl [( b{2) - Ci (2) ] = o . (89) 

From this, the solution for the total curvature perturbation at the first decay is 

Cl(2) = 4/ 70 i [C 7o (l) - &(!)] + 3/al [Ca(l) ~ Cl(l)] + 3/fcl [Cb(l) ~ Cl(l)] + /7olC 7o (2) + /alCo(2) + /fclCh(2) • (90) 



At second order, Eq. (I43|) reads 

« 7ll [4(c 7l(1) - Cm))] 2 + n 61 [3(c 6 (i) - Cm))] 2 + 4fi 71 i(C 7l( 2) - Ci( 2) ) + sn bl (c m - ( 1(2) ) = o, (91) 

which gives 

C 7l(2 ) = -4[C 71 (i)-Ci(i)] 2 + 3(l-i? 1 )[C b (i)-Ci(i)] 2 + Ci(2) + (l-i?i)C b (2)-(l-^i)Ci(2)- (92) 

2. Second decay 

We find the second order radiation perturbation after the second decay from Eq. (I44p [or from (|92p by first dropping 
curvaton b (which also implies i?i = 1), then relabeling a — > b, and 1 — > 2 and — ► 1], The result is 

C 72 (2) = -4 [C 72 (l) - C2(l)] 2 + C2(2) • (93) 

But recalling that £ 2 (i) = C 72 (i) this is just 

C2(2) = C 72 (2) , (94) 
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exactly as it should be since after the second decay only radiation is left. Now £2(2) is found from (190)) by first dropping 
curvaton b, and then relabeling a — > b, 1 — > 2 and — » 1. Thus we have 

C 2 (2) = 4(1 - / 62 ) [C 7l (l) ~ C2(l)] 2 + 3/ 62 [0,(1) - C2(l)] 2 + (1 - /b2)C 7 i(2) + /b2Cb(2) • ( 95 ) 

Here we substitute the curvature perturbations from our calculations above and assume 0y o (i) = Oy (2) = 0. Then we 
end up with the second order part of the primordial curvature perturbation 



where 



C(2) = C2(2) = C£ (1) + ^0 2 (1) + £0,(1)0(1) + FCa(2) + G0(2) , (96) 



C = - 2 -Rl/ol/62 - R lfalfb2 + 7R\ fa \fb2 ~ 4-Ri/ ol 

+ 3 -Rl/al ~ fal ~ Rlfal ~ Rlfblfal + 3i?l/al ~ 3/b2-Rl/ ol 

+fb2.fal + fblRlfal + fbzRlfblfal ~ 3/b2-Rl/al (97) 

D = -1 + f b 2Rifaifbi - 7Rihi - fbi + Vbi + 5Ri + f b 2 - ^Rlfbi + ZRffbi 
+7fb2Rifbi + 7R\fbifb2 - 14-Ri/bi/b2 - 2/ 2 2 i? 2 / 2 1 + 4/ 2 2 i? 2 /bi 
-fb2R\fbi + 2/ h 3 2 i?i/6i - 4J? X - 5/fc 2 i?i + 7/&2-R1 - 2/ 2 2 i? 2 - / b2 ^i 

+3-Ri/ fc 2 i - -Ri/bi + jte/bi - 2/ b2 / b i - Rifaifbi - 3/b2-Ri/fci + hiRifbi (98) 

■E = 2/qI + 2/b 2 i?l/al/bl - ZRlfalfbl ~ 10-Rl/al - ^Rlfblfal + 10/b2-Rl/al 
+ 2/b2-Rl/bl/ol + 14i?i/al/b2/bl — 4/ 2 2 i? 2 / a l/bl — 2/ 62 #l /al /bl 
— 6/b2-Rl/ol/bl + 8i? 2 / a i - 8Rlfalfbl - 14-R 2 /al/b2 + 4/ 2 2 i? 2 / a i 

+2fl 2 R{f al + QRlUhl - Valfbl - 2 /62/al + 2/62/al/bl (99) 

F = r a = {\-fm)f a xRx (100) 
G = r 6 = 1 - J?i + /bifli + / 62 -Ri - .fe/bi^i (101) 



In Appendix [51 we rewrite these coefficients with i?i from (158j) substituted in and organized in order of raising powers 
of /ai, /b2, and /bi. 

Assuming linear evolution between Hubble exit and the beginning of curvaton oscillation, the genuine second order 
curvaton perturbations in (|96|) read [see Eqs. ([17]) and (|18[) ] 

Ca(2) = -^Ca 2 ( l)' ( 102 ) 

0(2) = -fc? ( i). (103) 



Comparing now ([96f with (T2"l"j) we conclude 



C = C-|F, (104) 
D = D-\G. (105) 

(106) 

Recalling that A = r a and B = n from the first-order analysis, we now have found the non-linearity parameter /nl 
of Eq. ([H]): 

5 r 2 a (C - §r a ) + i/3 2 £r a rb + /3 4 rf (£> - f r 6 ) 

/nl = t; — 5 -. (107) 

6 (r 2 + /3 2 r 2 ) 2 



IV. /nl in various limits 



The general expression for /nl in Eq. (|107p written in terms of the coefficients C, D, E, r a and rb in Eqs. (|9"T)) - (|101|) 
[or equivalently ()B1|) (|B5[) in Appendix B] would be very complicated as a function of all four parameters f a i, fbi, 
fb2 and (3. Hence in this section we identify some useful limiting cases in which the expression significantly simplifies. 
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A. Single curvaton limits 

As a consistency check we first show that if either of the curvaton densities is negligible when they decay then the 
result (|107p simplifies to the single curvaton decay result [23|, H3, H3| 

fZ Kl V) = ^ - § - f , (108) 

where either / = f al if f b2 = 0, or / = f b2 if f a i = 0. 

If density of the second curvaton is always negligible, then we can set f b \ — fb2 = which leads to 

r a = fai, (109) 
C = 3/oi - 2/^ - / a 3 x , (110) 

L> = E = r b = 0, and then Eq. {TUT]) yields 

/NL = ^--^-^=/^ Sl °(/ a l). (Ill) 

On the other hand, if the density of the first curvaton is negligible, we have / al = which leads to 

n = fa, (112) 

D = 3/ b2 - 2/ b 2 2 - / b 3 2 , (113) 

C = E = r a = 0, and then 

/NL = ^-^-^ = /^ SlC (/ fc2 ). (114) 

Hence in the limit where only one curvaton has a non-negligible density, our two-curvaton result for /nl reduces to 
the well-known single-curvaton result (|108p . 



B. Simultaneous decay of curvatons a and b 

In the case where both curvatons decay at the same moment, we have f b \ = (1 + f a i/3)fb2, as discussed after 
Eq. (76]). Then we find 

r a = fai (115) 

r b = hi (116) 

C = -fal {fll + 2/„l + hlfal - 3) (117) 
D = -fhl (/hi + 2/fcl + / 6 l/al - 3) (118) 

£ = -2/„iAi(/ ol + / W +2). (119) 

The non-linearity parameter /nl which follows from these is indicated by the thin solid lines in Figs. [31 0] and [5] for 
the cases (3 — > cxd, /3 = 1, and /3 = respectively. 

In the limiting case where the second curvaton is homogeneous ((3 = 0) the expression (|107[) reduces to 

/NL = /^(/al) - jj/ M • (120) 

Conversely, in the limit when the first curvaton is homogeneous {(3 — > oo) we get 

/nl - fZ gl °(hi) - \hi • (121) 

Note that in either case the example of two curvatons which decay at the same time does not reduce exactly to the 
case of a single curvaton. However we can show that taking into account the constraint (|73p the minimum value for 
/nl (for any value of (3) is still that found for a single curvaton, min(/NL) = —5/4. 
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C. Both curvatons subdominant at decay 

If the energy density of both curvatons is small when they decay then f a \ <C 1 and fbi <C 1 [which implies fbi <C 1, 
from Eq. (176p ]. To first order in the / parameters, we find 



Ta fai (122) 

n fbi (123) 

C ~ 3/ ol (124) 

D ~ 3/ 62 (125) 

E ~ 0. (126) 



In this case the non-linearity parameter Eq. (| 1 0T|) reduces to the quasi-second order result (|88|) . As expected the 
quasi-second order result gives a good approximation when f a \ and /b 2 are small and the non-Gaussianity is thus 
large. 



D. Second curvaton negligible at the first decay 

If the energy density of the second curvaton b is negligible when the first curvaton a decays then we can set hi = 
and the linear result reduces to the simpler result (|T2[) found by Choi and Gong [l7|. We then have 

Ta = (I -fbi). fal (127) 

n = f b2 (128) 

C = hi (1 - hi) (3 - 2/„i - f al ) + f al hi (3 - 2f b i - Hi) (129) 

D = hi{Z-Vbi- fa) (130) 

E = -2/ a i(3/ fc2 -2/ fc 2 2 -/ fc 3 2 ) (131) 



and 

/NL = - / b2 ) 3 / a 3 l [| - 2/al ~ P al + hlfbl (hi + 3)] - [3 2 (1 - hlf fllftl (hi + 3) 

+/3 4 / b 3 2 (| - 2/ h2 _/&)}/{(!- / fe2 ) 2 / 2 x + /3 2 A 2 2 } 2 (132) 

In order to illustrate the general behavior of /nl we discuss three particular cases. 

a. Case (3 — > oo: this describes a situation where the first curvaton is effectively homogeneous, Ca = 0. This case 
is shown by the thick solid vertical lines in Fig. [3J 

In this case the universe remains homogeneous until the density of the second curvaton becomes non-negligible. 
Thus there is no difference from the original single inhomogeneous curvaton scenario, and Eq. (|132[) gives 

/nl = /^ sle (A 2 ) ■ (133) 

Note that we have assumed here that the the density of second curvaton is negligible when the first decays, hi = 0. 
If this is not the case then there can be a some dependence on f a \ (and hi) illustrated in Fig. [3] 

b. Case (3 = 1: this illustrates a situation where the fractional density perturbation in each curvaton field is 
comparable, here P$ a = P^ b . The thick solid lines in Fig. [4]show /nl(/&2, fai) corresponding to Eq. (|132|) with (3=1. 
We find that /nl can be large in this case only when fbi <C 1 and f a i <C 1 simultaneously. 

c. Case (3 = 0: describes a situation where the second curvaton is effectively homogeneous, (Cb = 0). With (3 = 
and hi = we have 



f 1 

JNL 



l-/( 



/,2 



/NL^C/al) + jj/M (Aa + 3) 



(134) 



In this case the inhomogeneous radiation produced by the first curvaton decay is diluted by the decay of the second 
homogeneous curvaton. The thick solid lines in Fig. [5] show fwhihi, fai) of Eq. (|134[) . 

When the second curvaton's density remains negligible throughout, hi <C 1, the non-linearity parameter is a 
function of a single variable f a i, and we recover the single-curvaton case (|111|) , which gives large /nl when f a i <C 1. 
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FIG. 3: /nl(//)2, fai) when f3 — * oo and f bl = (thick solid lines) or f bl = /t.2/2 (dotted lines) or /m = (1 + / a i/3)/&2 (t/i«n 
so^jrf lines). In the latter two cases the constraint (|73[) excludes a small region in the top right corner. Contours of equal /nl 
are shown, from left to right, for /nl = 1000, 114, 20, 5, 0, and — 1. 




FIG. 4: Same as Fig.[3]but now /nl(/m, fai) is shown for j3 = 1. Contours of equal /nl are shown, from the bottom left corner 
to the top right corner, for /nl = 1000, 114, 20, 5, 0, and — 1. 
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FIG. 5: Same as Fig. [3] but now fah{fa, fal) is shown for (3 = 0. Contours of equal /nl are shown, from bottom to top, for 
/nl = 1000, 114, 20, 5, 0, and -1. 



In the opposite limit where fbi ~* 1, we find a new regime in which the non-Gaussianity becomes large, which only 
appears in the presence of two (or more) decaying scalar fields. The non-linearity parameter is given by, 



1 



NL 



1-/62 



/nl S1c (/ Q i) 



10 

t 



(135) 



When f a i — > 1, the expression in square brackets approaches 25/12. Hence we have found an interesting result that 
/nl can be large even when / al ~ 1 if fa ~ 1. Note that the quasi-second order approximation (|88p correctly predicts 
/nl oc 1/(1 — fa) in this limit, but fails to reproduce the numerical coefficient 25/12. 



E. General case fa ^ 



As an example of the general case where the density of the second curvaton is non-negligible when the first curvaton 
decays (fa ^ 0), we show in Figs. [21 H] and [5] fah(fb2, fai) for the case where fa = fbi/^ and for the simultaneous 
decay where fbi = (1 + f a i/3)fa- Generally we find only a weak dependence on fbi-, although the constraint 
fbi < 1 — fai excludes a small region in the top right corner where fbi ~ 1 and f a i ~ 1. For intermediate values of 
fbi 6 [0, (1 + / a i/3)/b2] we find that the values of /nl fall between the cases fbi = and fbi = (1 + / a i/3)/b2, i-e., 
the contours of equal values of /nl in Figs. [3l [4] and [5] lie between the thick and thin solid lines shown. 

Finally, we find numerically that the minimum value of the nonlincarity parameter (|10T|) is min(/NL) = — 5/4, 
when the constraints (|73[) and (|76[) are taken into account. This is the same as the minimum value found for a 
single curvaton, min(/^' L Lslc ), assuming that the evolution between Hubble exit and beginning of curvaton oscillation 
is linear, which is the case, for example, for a weakly interacting field with a quadratic potential. 



V. CONCLUSION 



In this paper we have calculated the non-linear primordial curvature perturbation, following the decay of two 
curvaton fields in the early universe. The full non-linear curvature perturbation is given by the series of expressions, 
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(|4Tj) . (|43p. (J44J) and (|45|) . relating the final perturbation in the radiation after the second curvaton decays, £ 72 , to the 
initial curvature perturbations in each curvaton field, £ a and in the sudden-decay approximation. Expanding to 
first- and then second-order in the perturbations we have derived the general expressions in Eq. (|71| for the power 
spectrum and Eq. (|107p for the non-linearity parameter, /nl, which describes the non-Gaussianity of the primordial 
curvature perturbation at leading order. 

Large values of the non-linearity parameter in the two-curvaton model can qualitatively be understood by a sim- 
plified quasi-second-order model 



c 



Co(l) 



[Ca(l) 



Cb(l 



3,2 

4^(1) 



(136) 



leading to 



f: 



NL 



5 rl+l3 A i 



(137) 



where r a and Tb are given by first-order results (|66[) and 

We easily recover single field results: /nl ~ 5/(4r(,) if f a i = 0, or /nl ~ 5/(4r a ) if fbi = and fbi = 0. In either 
case the non-linearity becomes large if the transfer of curvaton to curvature perturbation is inefficient (r a <C 1 and 
rb -C 1). More generally in the two-curvaton model we find large /nl if both curvatons are subdominant at their 
decay time. 

We can also obtain large non-Gaussianity in a different way, only possible due to the existence of two curvatons, 
where the second curvaton is effectively homogeneous: (3 = 0. In this case /nl <x l/r a , but the efficiency r a oc 
/ a i(l — fb2) becomes small either when f a i is small (as in the one curvaton case) or when fb2 ~ 1- In this latter case 
the inhomogeneous radiation produced by the first curvaton decay is diluted by the decay of the second homogeneous 
curvaton. 

In all cases we find /nl > — 5/4, which seems to be a robust lower bound in single and multi-field curvaton models in 
which the curvaton field perturbations are themselves Gaussian, which should be a good approximation for a weakly 
interacting field. 
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APPENDIX A: NOTATION 



The notation becomes quite complicated with the need to specify two different curvaton species and the radiation 
produced in their decays, which happen at two different moments. Furthermore, we need to keep track on the first 
and second order parts of perturbations. Therefore, to advice the reader, we present in this section a list of some 
symbols that appear in our calculations and results. 

A subscript 1 or 2 after specification of the particle species tells whether the quantity is evaluated at first or sec- 
ond decay. Note that some quantities are constant, and hence don't need this subscript as apparent from the list below. 

The energy density parameters (fli = Pi/ ptot&i)'- 

fl a i the density of the first curvaton a just before the first decay, 

the density of the second curvaton b at the first decay, 
f2f,2 the density of the second curvaton b just before the second decay, 
51 7o i the density of pre-existing radiation at the first decay, 
fi 7l i the density of all radiation immediately after the first decay, 
f2 7l 2 the density of all radiation just before the second decay, 
f2 72 2 the density of all radiation immediately after the second decay. 

The full non-linear curvature perturbations: 

C (— C2) the primordial perturbation after the second curvaton decay, but before the nucleosynthesis, 
Ci the total perturbation at the first decay, 
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(2 the total perturbation at/after the second decay, 

the perturbation of the first curvaton a, 
(b the perturbation of the second curvaton b, 
C 7o the pre-existing radiation perturbation, 
£ 7l the radiation perturbation after the first decay, 
C72 ( = C2 = C) the radiation perturbation after the second decay. 

Subscript in parenthesis after all the other subscripts denotes the order considered. For example: 
£m the first order part of the primordial perturbation, 
£(2) the second order part of the primordial perturbation, 
Cim the first order part of the total perturbation at the first decay, 
£1(2) the second order part of the total perturbation at the first decay, 
£2(1) the first order part of the total perturbation at the second decay, 
£2(2) the second order part of the total perturbation at the second decay, 
Call) the first order part of the first-curvaton perturbation, 
Ca(2) the second order part of the first-curvaton perturbation, 
Cb(i) the first order part of the second-curvaton perturbation, 
Cb(2) the second order part of the second-curvaton perturbation, 
C70 
C70 
C71 
C71 
C72 
C72 



the first order part of the pre-existing radiation perturbation, 
2) the second order part of the pre-existing radiation perturbation, 
!) the first order part of the radiation perturbation after the first decay, 
2) the second order part of the radiation perturbation after the first decay, 

1) the first order part of the radiation perturbation after the second decay, 

2) the second order part of the radiation perturbation after the second decay. 



The energy density ratios (curvature perturbation transfer efficiencies) at the first decay: 

f _ 3Sl al 

■> al 4^, i+3O al +3O i)1 ' 
f _ 3^bi 

Jbl 4n T0 i+3O a i+3n 6 i ' 
f 1 

J7ol 4Q T0 i+3£lai+3ttbi ' 
r> _ 4-Sl b i 
1X1 4-4f2 61 ' 

The energy density ratios (curvature perturbation transfer efficiencies) at the second decay: 

3n b2 



■f b2 4n 7l2 +30 6 2 ' 

p 4£7 71 2 

^Ti 2 ~~ 40 T , 2 +3f2 b2 ' 



The total first order curvaton perturbation transfer efficiencies: 

„ _ (l-/b2)(3+/ a l)/„l 

1(1 3(l-/ bl )+/ a i ' 

„ _ (l-/6l)/b2(3+/al) + /6l/al 
' b ~ 3(1-/5! ) + 

The non-linearity parameters: 

the non-linearity parameter of the first curvaton perturbation £ a , 
the non-linearity parameter of the second curvaton perturbation 

/nl the non-linearity parameter of the primordial perturbation £ = £2 = C72 • 



APPENDIX B: COEFFICIENTS IN THE SECOND ORDER PART OF C 

In this section we express the coefficients C, D, E, F and G, Eqs. (|97til01[> . with Ri from Eq. (|58|) substituted 
in. We don't write the neatest possible forms, but instead organize the results in order of raising powers of of f a i, 
fb2, and fbi- This is particularly practical for the purposes of Sect. ITV1 where we express /nl in various limits. For 
example in the limit where both curvaton are subdominant at their decay time, i.e., f a i 1, fb2 <C 1, and fbi <C 1. 



The coefficients in the second order part of £ in (|96|) are: 
1 



C 



D 



E 



F 
G 



{ [27 - 27/ 61 



(3 - 3/ 61 + / Ql ) 2 

27/ 61 + (27 + 27/ bl )/ b2 

f2 , foe of2 



18/ 6 2 2 ~ 9/ 3 2 
18 + 3/& + (36 - 3/&)/w - 12/ b 2 2 - 6/ b2 



(-27 + 27/ 6 i)/ w j/„i 

/al 

/al 



+ 



8 + 2/ H + (ll-2/ 61 )/ ( 



62 



2/6 2 2 



/&2 



/al 



1 + /62 
1 



/ a 5 r} 



(3 - 3/ M + /al) 2 

+ 9/bi - + oj bl -f 1,10 - 4o/ fc i 
+ (-6 + 12/ bl -6/ 2 1 )/ 3 J/ai 



{(27 - 54/6i + 27/ 2 1 )/ b2 + (-18 + 36/ M - ISftJfa + (-9 + 18/ 6 i - 9/&)/, 



3/4 + (18 - 45/w + 30/& - 3/&)/w + (-12 + 24/ bl - 12/&)/? 



62 



3/61 - 8/ b 2 ! + 2/ b 3 x + (3 - 9/ M + ll/ 2 ! - 2/30/62 + (-2 + 4/ M - 2/ 2 1 )/ 2 2 + (-1 + 2/ M - 



A 2 ! + /62/&1 



(3 - 3/fti + /al) 2 



/al} i 



I8/61 - 18/ 



-54 + 36/6i + 18/6^/62 + (36 - 36/ fc i)/ b 2 2 + (18 - 18/ bl )/ b 3 2 



24/61 + 6/ 3 ! + (-36 + 6O/61 - 6/^/62 + (24 - 24/ b i)/ b 2 2 + (12 - 12/ bl )/ 3 2 f al 



+ 
+ 

r a 

n 



16/61 + 4/& + (-6 + 22/ bl - 4/ 2 1 )/62 + (4 - 4/6i)/ 2 2 + (2 - 2/ bl )/ 3 2 



f 3 

Jal 



2/&i + 2/bi/i 



62 



/a\} 



, f S >< {(3 - 3/6 2 )/al + (1 - /62)/a 2 l| , 
(o - 3/fti + Jal) } 

7^ T , v X {(3 - 3/ bl )/ b2 + [/61 + (1 - fa) fa] fal} ■ 

(3 — 3/61 + Jal) <■ > 
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